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We present an amplitude analysis of the decay D? > K-ntntn- based on a data sample of 


2.93 fb^! acquired by the BESIII detector at the (3770) resonance. With a nearly background 
free sample of about 16000 events, we investigate the substructure of the decay and determine the 
relative fractions and the phases among the different intermediate processes. Our amplitude model 
includes the two-body decays D? > K*°p®, D? => K- af (1260) and D? > Kj (1270)z*, the three- 
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body decays D? > K*°nta7~ and D? > Korto, 


D?  K^n*m«m 
fraction of 54.696. 


PACS numbers: 13.20.Ft, 14.40.Lb 
I. INTRODUCTION 


The decay D? > K-ntata- is one of the three 
golden decay modes of the neutral D meson (the other 
two are D? > K^«* and D?  K^«*«?). Due to a 
large branching fraction and low background it is well 
suited to use as a reference channel for other decays of 
the D? meson [1]. An accurate knowledge of its resonant 
substructure and the relative amplitudes and phases are 
important to reduce systematic uncertainties in analy- 
ses that use this channel for reference. In particular, 
the lack of knowledge of the substructure leads to one of 
the largest systematic uncertainties in the measurement 
of the absolute branching fractions of the D hadronic 
decays [. The knowledge of the decay substructure in 
combination with a precise measurement of strong phases 
can also help to improve the measurement of the CKM 
angle y (the phase of Va relative to Vu») [3. In the 
measurement of 7, the parametrization model is an im- 
portant input information in a model dependent method 
and also can be used to generate Monte Carlo (MC) 
simulations to check the sensitivity in a model indepen- 
dent method 4. Furthermore, the branching fractions 
of intermediate processes can be used to understand the 
D? — D? mixing in theory 


The decay D? > K-«z*«*-- was studied by Mark 
III M and E691 [8] more than twenty years ago. Both 
measurements are affected by low statistics. Using about 
1300 signal events, Mark III obtained the branching frac- 
tions for D" => Kat (1260), D" — K*955 Dp? => 
K;(1270)r*, as well as for the three- and four-body 
nonresonant decays. Based on 1745 signal events and 
800 background events, E691 obtained a similar result 
but without considering the D? > KI (1270)m+ decay 
mode. The results from Mark III and E691 have large 
uncertainties. Therefore, further experimental study of 
D? + K-«*«*n- decay is of great importance for im- 
proving the precision of future measurements. 


In this paper, a data sample of about 2.93 fb ^! [g] 
collected at the (3770) resonance with the BESIII detec- 
tor in 2010 and 2011 is used. We perform an amplitude 
analysis of the decay D? — K^ x*«*«- (the inclusion of 
charge conjugate reactions is implied) to study the reso- 
nant substructure in this decay. The (3770) decays into 
a D? D? pair without any additional hadrons. We employ 
a double-tag method to measure the branching fraction. 
In order to suppress the backgrounds from other charmed 
meson decays and continuum (QED and qq) processes, 
only the decay mode D? > K*«- is used to tag the 
D? D? pair. A detailed discussion of background can be 
found in Sec. The amplitude model is constructed 


as well as the four-body nonresonant decay 


-. The dominant intermediate process is D^ — K~ aj (1260), accounting for a fit 


using the covariant tensor formalism [11]. 


II. DETECTION AND DATA SETS 


The BESIII detector is described in detail in Ref. fg. 
The geometrical acceptance of the BESIII detector is 
93% of the full solid angle. Starting from the interaction 
point (IP), it consists of a main drift chamber (MDC), a 
time-of-flight (TOF) system, a CsI(T1) electromagnetic 
calorimeter (EMC) and a muon system (MUC) with lay- 
ers of resistive plate chambers (RPC) in the iron return 
yoke of a 1.0 T superconducting solenoid. The momen- 
tum resolution for charged tracks in the MDC is 0.5% at 
a transverse momentum of 1 GeV/c. 


Monte Carlo (MC) simulations are based on 
GEANTA [13]. The production of (3770) is simulated 
with the KKMC package, taking into account the 
beam energy spread and initial-state radiation (ISR). 
The PHOTOS package is used to simulate the final- 
state radiation (FSR) of charged tracks. The MC sam- 
ples, which consist of (3770) decays to DD, non-DD, 
ISR production of low mass charmonium states and con- 
tinuum processes, are referred to as “generic MC” sam- 
ples. The EvtGen package is used to simulate the 
known decay modes with branching fractions taken from 
the Particle Data Group (PDG) ini and the remain- 
ing unknown decays are generated with the LundCharm 
model [17. The effective luminosities of the generic MC 
samples correspond to at least 5 times the data sample 
luminosity. They are used to investigate possible back- 
grounds. The decay D? > K(r*w )K 7^ has the 
same final state as signal and is investigated using a ded- 
icated MC sample with the decay chain of (3770) > 
D°D° with D} > RR EU and D? + K*m-, re- 
ferred to as the “K2Ka MC". The decay model of 
D? => K3K~x* is generated according to CLEO’s re- 
sults [18]. In amplitude analysis, two sets of signal 
MC samples using different decay models are generated. 
One sample is generated with an uniform distribution 
in phase space for the D? > K-n+a*n- decay, which 
is used to calculate the MC integrations and called the 
*PHSP MC” sample. The other sample is generated ac- 
cording to the results obtained in this analysis for the 
D? => K-ntatan- decay. It is used to check the fit 
performance, calculate the goodness of fit and estimate 
the detector efficiency, and is called the *SIGNAL MC” 
sample. 
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III. EVENT SELECTION 


Good charged tracks are required to have a point of 
closest approach to the interaction point (IP) within 10 
cm along the beam axis and within 1 cm in the plane 
perpendicular to beam. The polar angle 0 between the 
track and the e* beam direction is required to satisfy 
[cos0| < 0.93. Charged particle identification (PID) 
is implemented by combining the energy loss (dE/dx) 
in the MDC and the time-of-fight information from the 
TOF. Probabilities P(K) and P(r) with the hypothe- 
ses of K or « are then calculated. Tracks without PID 
information are rejected. Charged kaon candidates are 
required to have P(K) > P(r), while the m candidates 
are required to have P(t) > P(K). The average efficien- 
cies for the kaon and pions in K z^ 7*7 are ~ 98% and 
~ 99% respectively. The D9D9 pair with D? > Ktn- 
and D? > K ^a *«^m- is reconstructed with the require- 
ment that the two D? mesons have opposite charm and 
do not have any tracks in common. Since the tracks in 
K «'m«m- have distinct momenta from those in Km, 
misreconstructed signal events and K/7 particle misiden- 
tification are negligible. Furthermore, a vertex fit with 
the hypothesis that all tracks originate from the IP is 
performed, and the x? of the fit is required to be less 
than 200. 


For the Ktm and K^«'*-«*- combinations, two 
variables, Mpc and AE, are calculated: 


MBc = V El cam iar Pp: (1) 


and 
AE= Ep ii beam (2) 


where pp and Ep are the reconstructed momentum and 
energy of a D candidate, Epeam is the calibrated beam 
energy. The signal events form a peak around zero in 
the AE distribution and around the DÓ mass in the 
Mnpo distribution. We require —0.03 < AE < 0.03 GeV 
for the K *«- final state, —0.033 < AE < 0.033 GeV 
for the K~atata~ final state and 1.8575 < Mpa < 
1.8775 GeV/c? for both of them. The corresponding AE 
and Mgc of selected candidate are shown in Fig.[I] where 
the background is negligible. 


To ensure the D? meson is on shell and improve the 
resolution, the selected candidate events are further sub- 
jected to a five-constraint (5C) kinematic fit, which con- 
strains the total four-momentum of all final state parti- 
cles to the initial four-momentum of the ete" system, 
and the invariant mass of signal side K7~atata7 con- 
strains to the D? mass in PDG [I]. We discard events 
with a x? of the 5C kinematic fit larger than 40. In or- 
der to suppress the background of D? > KOK s^ with 
KÌ => ata, which has the same final state as our signal 
decay, we perform a vertex constrained fit on any 7? 77 
pair in the signal side if the ++ c^ invariant mass falls into 


the mass window |m,+,- —mgo| < 0.03 GeV/c? (mio is 
the K? nominal mass H), and reject the event if the cor- 
responding significance of decay length (e.g. the distance 
of the decay vertex to IP) is larger than 20. The K$ veto 
eliminates about 80% D? + K9K r* background while 
retaining about 9996 of signal events. After applying all 
selection criteria, 15912 candidate events are obtained 
with a purity of 99.496, as estimated by MC simulation. 


The MC studies indicate that the dominant back- 
ground arises from the D? + K9K «* decay, the cor- 
responding produced number of events is estimated ac- 
cording to 


N(KSK -v'|K*m )— 
Y(K^m'n'm |K*m P. 
e(K-m*mtmn-|K*n-) 


B(K9K-m*) (3) 
B(K-m*mw*tmn-) 


where N(KgK m-^*|K*m- ) is the production of 
v(3770 => DD? with D? > KQK--*' and 
D? — Ktm, Y(K-n*n*m-|K*m-) is the signal 
yield with background subtracted but without efficiency 
correction applied and e is the corresponding efficiency 
obtained from the SIGNAL MC sample, which is gener- 
ated according to the results of fit to data whose peaking 
background estimated from the generic MC sample. 
B(K-ntatn—-) and B(KgK m^) are the branching 
fractions for D? > K-atata- and D? > K9K-s*, 
respectively, which are quoted from the PDG [1]. Ac- 
cording to Eq. (3), the number of peaking background 
events (Npeaking) is estimated to be 96.8 + 14.5. 


All other backgrounds from DD, qq and non-DD de- 
cays are studied with the generic MC sample. Their to- 
tal contribution is estimated to be less than ten events, 
of which 5.5 and 2.0 are from the D?D? decays and 
the non-DD decays, respectively. These backgrounds 
are neglected in the following analysis and their effect 
is considered as a systematic uncertainty, as discussed in 


Sec. 


IV. AMPLITUDE ANALYSIS 


The decay modes which may contribute to the D? > 
K-n«*m-*m-- decay are listed in Table[I] where the sym- 
bols S, P, V, A, and T denote a scalar, pseudoscalar, 
vector, axial-vector, and tensor state, respectively. The 
letters S, P, and D in square brackets refer to the relative 
angular momentum between the daughter particles. The 
amplitudes and the relative phases between the different 
decay modes are determined with a maximum likelihood 
fit. 
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FIG. 1. Distributions of data for AE [(a) and (c)] and Mao [(b) and (d)] in K*-^ side [(a) and (b)] and in K^ «*^-*«- side 
[(c) and (d)]. The arrows indicate the selection criteria. In each plot, all selection criteria described in this section have been 


applied except the one on the variable. 


A.  Likelihood function construction 


The likelihood function is the product of the probabil- 
ity density function (PDF) of the observed events. The 
signal PDF fs(p;) is given by 


e(pi)|M (p;)| Ra(p;) 
J e(pi)|M (pj) |? Ra(p;)dp; ' 


where e(p;) is the detection efficiency parametrized in 
terms of the final four-momenta p;. The index j refers 
to the different particles in the final state. R4(p;)dp; is 
the standard element of the four-body phase space 
which is given by 


fs(p;) = (4) 


9 


4 
Ra(p;)dp; = ôt | ppo — 
J 


4 
si- ©) 
? | 11 (2x)?25E; 
1 j—1 


M(p;) is the total decay amplitude which is modeled as 
a coherent sum over all contributing amplitudes 


M(p;) =Y cn As(pi), (6) 


where the complex coefficient c, = py eio» (pn and dy are 
the magnitude and phase for the n*" amplitude, respec- 
tively) and An (p;) describe the relative contribution and 
the dynamics of the n*^ amplitude. In four-body decays, 
the intermediate amplitude can be a quasi-two-body de- 
cay or a cascade decay amplitude, and A,(p;) is given 
by 


An(pj) = Pi (mi) P? (ma)8s (pj) Fi (pj) F2 (pj) FP (n5).(7) 


where the indices 1 and 2 correspond to the two interme- 
diate resonances. Here, P?(m4) and F?(p;) (a = 1,2) 


are the propagator and the Blatt-Weisskopf barrier fac- 
tor [19], respectively, and F? (p;) is the Blatt- Weisskopf 
barrier factor of the D? decay. The parameters m; and 
mg in the propagators are the invariant masses of the 
corresponding systems. For nonresonant states with or- 
bital angular momentum between the daughters, we set 
the propagator to unity, which can be regarded as a very 
broad resonance. The spin factor S,(p;) is constructed 
with the covariant tensor formalism . In practice, the 
presence of the two 7* mesons imposes a Bose symmetry 
in the K^«*7-*-«- final state. This symmetry is explic- 
itly accounted for in the amplitude by exchange of the 
two pions with the same charge. 


'The contribution from the background is subtracted in 
the likelihood calculation by assigning a negative weight 
to the background events 


Naata Nokg 
In L = 5 In fs(p7) — 5 wy In fs(p} J; (8) 
k=1 k'=1 


where Naata is the number of candidate events in data, 
were and Nyyg are the weight and the number of 
events from the background MC sample, respectively. 
In the nominal fit, only the peaking background D? => 
KO K-z* is considered, and the weight wp is fixed to 
Npeaking/Nokg- pr and pF are the four-momenta of the 
jt} final particle in the k*" event of the data sample and 
in the k*t! event of the background MC sample, respec- 
tively. 


The normalization integral is determined by a MC 
technique taking into account the difference of detector 
efficiencies for PID and tracking between data and MC 
simulation. The weight for a given MC event is defined 
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as 


€j.data(Dj) 
*« (pj) II €;.Mc(D;) , (9) 
where €jdata(pj) and ej c(p;) are the PID or track- 
ing efficiencies for charged tracks as a function of p; for 
the data and MC sample, respectively. The efficiencies 
€j,data(pj;) and ej,uc(pj) are determined by studying the 
D? + K~xtrtr sample for data and the MC sample 
respectively. The MC integration is then given by 


/ e(p;)|M (pj)? Ra(p;)dp; 


1 
7 1 Nwmc |M (p~s )) 24, (p~s) ( 0) 
Nc 2 | Mee (ph yt 


kmc 


, 


where kmc is the index of the Es event of the MC sam- 


ple and Nmc is the number of the selected MC events. 
M?" (p;) is the PDF function used to generate the MC 
samples in MC integration. In the numerator of Eq. (4), 
e(p;) is independent of the fitted variables, so it is re- 
garded as a constant term in the fit. 


1. Spin factors 


Due to the limited phase space available in the decay, 
we only consider the states with angular momenta up to 
2. As discussed in Ref. [11], we define the spin projection 


E 
operator [o NM for a process a — bc as 


PO) = —guv» + Tara (11) 
for spin 1, 
2 = 
P unt = 
1 (12) 


Loa) pt) ü np. Lp) p) 
SD AT a e PLUR 


1417 H2V2 Hiva' bevy 1H27 viva 


for spin 2. The covariant tensors is ies for the final 
states of pure orbital angular momentum L are con- 
structed from relevant momenta pa, Pb, Pe 


rino = (D Peat, (18) 
where r = pp — pe. 


Ten kinds of decay modes used in the analysis are listed 
in Table [I] We use TU as to represent the decay from 


the D meson and ïo), pz to represent the decay from the 
intermediate state. 


2. Blatt-Weisskopf barrier factors 


The Blatt-Weisskopf barrier factor Fr(pj) is a 
function of the angular momentum L and the four- 
momenta pj; of the daughter particles. For a process 


à — bc, the magnitude of the momentum q of the daugh- 
ter b or c in the rest system of a is given by 


Sa + Sb - 8c)? 
j= es ; (14) 


with sg = E2 — ph, = a,b,c. The Blatt-Weisskopf 
barrier factor is then given by 


Fr(q) = z^ Xy (q), (15) 


where z — qR. R is the effective radius of the barrier, 
which is fixed to 3.0 GeV! for intermediate resonances 
and 5.0 GeV ^! for the D? meson. Xz (q) is given by 


Xr-o(q) = 1, (16) 


X1=1(9) = 1/ E (17) 
Xil) = zr (2) 


3. Propagator 


The resonances K*? and aï (1260) are parametrized as 
relativistic Breit- Wigner function with a mass depended 
width 


1 


Pn (m2 = sa) — imoL'(m)* 


(19) 


where mo is the mass of resonance to be determined. 
T(m) is given by 


2L41 2 
X 
rem =m (Z) (Œ) (FEF). e 
qo m / \ Xr(qo) 
where qo denotes the value of q at m = mo. The 
Ky (1270) is parametrized as a relativistic Breit-Wigner 
function with a constant width T(m) = To, and the p? is 


parametrized with the Gounaris-Sakurai line shape [20], 
which is given by 


r rm A" at 
Pas (m) (m2 — m2) + f (m) = imor (m) (21) 


where 


f(m) = To 2 [g2(h(m) = h(mo)) 
do 


T (22) 
2 2) 2 
TURO m )4 d(m?) "TR 
and the function h(m) is defined as 
2 2 
h(m) = Í m (= 2), (23) 
T m 2Mr 
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TABLE I. Spin factors S(p) for different decay modes. 


Decay mode S(p) 
^ DIS|— V,V>5.V, S E.P». Vo S PP. O Aan (NS) ` C;C;OC*” 
D|S] > Vi V2, Vi > Pı P2, Vo > P3 P4 COP Vi JEP (Va) 


D|P] — Vi V2, Vi => Pi Pa, Vo => P3P4 
D[D] > V4 Va, Vı > Pi P3; V2 +> P3P4 
D > APi, A[S] > VP5, V > P3P4 

D => AP), A[D| > VPs, V — Pa Pa 

D => AP1, A > SP2, S 2 Ps P4 

D > VS, V => P, Po, S > Ps P4 

D ViPi, Vi V2 P2, V2 > P3 P4 

D > PPi, Po VP», V => P3Pa 

D > TS, T Pi Pa, S > P3P4 


with 
dh 


2 2 1] 
d(m2) m2=më Í mal 


= h(mo)[(843) `` — 
where My is the charged pion mass. The normal- 


ization condition at Pas(0) fixes the parameter d = 
f(0)/ (Tomo). It is found to be 


3 m2 " (= zm) 4 mo m2mo 
27 qo T qÈ 


4. Parametrization of the Kr S-wave 


For the Kr S-wave [denoted as (K7)g wave], we use 
the same parametrization as BABAR 5 which is ex- 
tracted from scattering data La). The model is built from 
a Breit-Wigner shape for the K3(1430)? combined with 
an effective range parametrization for the nonresonant 
component given by 


A(MKr+) F sinópe??F + Rsin dpe? e'?F , (26) 
with 
óp = or cot"! |Z + TE (27) 
ee aq 2 |? 
— 1 T 
ÒR = Ọr + tan a l (28) 


where a and r denote the scattering length and effec- 
tive interaction length. F (ór) and R (ón) are the rel- 
ative magnitudes (phases) for the nonresonant and res- 
onant terms, respectively. q and I (mgr) are defined as 
in Eq. (I4) and Eq. (20), respectively. In the fit, the 
parameters M, T, F, óp, R, or, a and r are fixed to 
the values (Lined from the fit to the T K8ntn 
Dalitz plot [21], as summarized in Table [I | These fixed 
parameters eri be varied within their uncertainties to es- 
timate the corresponding systematic uncertainties, which 
is discussed in detail in Sec. 


+ (27m?) ) (24) 


pis 


Envie p" (DITO (DEO V E®? (V3) 
UM dA a V3) 
TED) Pi (A)EO”(V) 
POUDES (A) 

TOM DP (V) 

m v c 
E d 
al (Pa)l (V) 
TOR (DA ( T) 


TABLE II. Kr S-wave parameters, obtained from the fit to 
the D? > KOz*z- Dalitz plot from BABAR [21]. 
M(GeV/c?) 1.463 + 0.002 


T(GeV/c?) 0.233 + 0.005 


F 0.80 + 0.09 
Ór 2.33 + 0.13 
R 1 (fixed) 
or —5.31 + 0.04 
a 1.07 + 0.11 
r —1.8 + 0.3 


B. Fit fraction and the statistical uncertainty 


We divide the fit model into several components ac- 
cording to the intermediate resonances, which can be 
found in Sec. The fit fractions of the individual com- 
ponents (amplitudes) are calculated according to the fit 
results and are compared to other measurements. In 
the calculation, a large phase space (PHSP) MC sample 
with neither detector acceptance nor resolution involved 
is used. The fit fraction for an amplitude or a component 
(a certain subset of amplitudes) is defined as 


er [An (o) 
FF(n) = E» 29 
m Xx MP? ad 


where Ån(p¥) is either the n'" amplitude [An(p¥) = 
c» An(p*)] or the nt^ 
amplitudes [An(p*) = V ith; An; (p7)], Ngen is the num- 
ber of the PHSP MC events. 


component of a coherent sum of 


To estimate the statistical uncertainties of the fit frac- 
tions, we repeat the calculation of fit fractions by ran- 
domly varying the fitted parameters according to the er- 
ror matrix. Then, for every amplitude or component, we 
fit the resulting distribution with a Gaussian function, 
whose width gives the corresponding statistical uncer- 
tainty. 
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C. Goodness of fit 


To examine the performance of the fit process, the 
goodness of fit is defined as follows. Since the D? and 
all four final states particles have spin zero, the phase 
space of the decay D? > Kortnr can be completely 
described by five linearly independent Lorentz invariant 
variables. Denoting as TI the one of the two identical 
pions which results in a higher rtr invariant mass and 
the other pion as 77 , we choose the five invariant masses 
m + e Tuo go gt ao Mr tate _ and Te + at: To 
calculate the goodness of fit, the five- dimensional: phase 
space is first divided into cells with equal size. Then, ad- 
jacent cells are combined until the number of events in 
each cell is larger than 20. The deviation of the fit in each 
Np— 
quantified as x? = 375 4 xl, where N, and NS*P are the 
number of the observed events and the expected number 
determined from the fit results in the p/^ cell, respec- 
tively, and n is the total number of cells. The number of 
degrees of freedom (NDF) v is given by v = (n—1)—npa:, 
where npa: is the number of the free parameters in the 
fit. 


cell is calculated, xp = and the goodness of fit is 


V. RESULTS 


In order to determine the optimal set of amplitude that 
contribute to the decay D? > K-^«*«*-- , considering 
the results in PDG [1], we start with the fit including the 
components with significant contribution and add more 
amplitude in the fit one by one. The corresponding sta- 
tistical significance for the new amplitude is calculated 
with the change of the log-likelihood value Aln L, taking 
the change of the degrees of freedom Av into account. 


In the Krt and ntr” invariant mass spectra, there 
are clear structures for K*° and p°. The intermediate res- 
onance K, (1270) is observed with K7 (1270) > K*?z— 
or K-p?. In the ttata- invariant mass spectrum, a 
broad bump appears. We find this bump can be fitted as 
aï (1260), which was also observed by the Mark IIT [ai ex- 
periment. If it is fitted with a nonresonant (p°7*), am- 
plitude instead, we find that the significance for aï (1260) 
with respect to (p?m^)A is larger than 10c. The three- 
body nonresonant states come from two kinds of con- 
tributions, Krt? and K? ntr. The K*95-/ K - o? 
can be in a pseudoscalar, a vector or an axial-vector 
state, while the K~a*/n*a~ can be in a scalar state. 
The four-body nonresonant states are relatively complex, 
such as D > VV, D => VS, D > TS, D > TV, D > 
AP with A — VP or SP, all of which may contribute 
to the decay. Since the process D? => K-a[(1260), 
aï (1260)[S] > p°xt has the largest fit fraction, we fix 
the corresponding magnitude and phase to 1.0 and 0.0 
and allow the magnitudes and phases of the other pro- 


cesses to vary in the fit. 


We keep the processes with significance larger than 
5o for the next iteration. The fit involving both the 
Kat (1260) and the nonresonant K~(p°x*), contri- 
bution does not result in a significantly improvement 
of fit, but the fit fractions of the two amplitudes are 
much different with the assumption of only K~ af (1260) 
and are nearly 100% correlated. We avoid this kind 
of case and only consider the resonant term, in agree- 
ment with the analysis of Mark III Gi For the process 
D? => Kj (1270)r+ with Kj (1270)[S] > K*?z-, the 
corresponding significance is found to be 4.30 only, but 
we still include it in the fit since the corresponding D- 
wave process is found to have a statistical significance of 
larger than 9c. Better projections in the invariant mass 
spectra and an improved fit quality x? are also seen with 
this S-wave process included. 


Finally, we retain 23 processes categorized into seven 
components. The other processes, not used in our nom- 
inal results but have been tested when determining the 
nominal fit model, are listed in Appendix A. The widths 
and masses of K*? and p? are determined by the fit. The 
results of are listed in Table The Kj (1270) has a 


TABLE III. Masses and widths of intermediate resonances 
K*? and p°, the first and second uncertainties are statistical 
and systematic, respectively. 


Resonances Mass (MeV/c?) Width (MeV/c?) 
K* 894.78 £0.75 + 1.66 44.18 + 1.57 + 1.39 
p? 779.14 + 1.68 + 3.98 148.42 + 2.87 + 3.36 


small fit fraction, and we fix its mass and width to the 
PDG values [1]. The af (1260) has a mass close to the up- 
per boundary of the c *7*«- invariant mass spectrum. 
'Therefore, we determine its mass and width with a like- 
lihood scan, as shown in Fig.[2] The scan results are 


M+ (1260) = 1362 +13 MeV/c’, 


À (30) 
T a+ (1260) =542 +29 MeV/e?, 


where the uncertainties are statistical only. The mass 
and width of af (1260) are fixed to the scanned values in 
the nominal fit. 


Our nominal fit yields a goodness of fit value of x? /v = 
843.445/ 748 = 1.128. To calculate the statistical signifi- 
cance of a process, we repeat the fit process without the 
corresponding process included, and the changes of log- 
likelihood value and the number of free degree are taken 
into consideration. 'The projections for eight invariant 
mass and the distribution of x are shown in Fig.[3] All 
of the components, amplitudes and the significance of 
amplitudes are listed in Table 'The fit fractions of 
all components are given in Table[V] The phases and fit 
fractions of all amplitudes are given in Table [VI] 
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FIG. 2. Likelihood scans of the width (a) and mass (b) of aj (1260). 
TABLE IV. Statistical significances for different amplitudes. 
Component Amplitude Significance (co) 
- DS] > Kp > 10.0 
D? + K*959 D?[P] ^ Ro, 2 > 10.0 
D?[D| > K*°p K > 10.0 
5 D" > K- a1 (1260), a} (1260)[S] > p°x > 10.0 
0 + + D „+ 1 , d 
Herde a ON au (1260) cam D? = AE t(1260)[D] > p*r + 7.4 
= D? > Kj (1270) F, Kj (1270)[S] > Kx 4.3 
0 pam au = x0, — 1 = A 
D > Kï UU. Kï (1270) > K" 77 po, => (1270), Kj (1270)(D ] > Rn 9.6 
D" => K; (1270)^, Ky (1270) > K^ p D" => Ky (1270)^, Kj (1270)|$] > K^ p > 10.0 
DY > (pK )am TOR )a[D] > K^ p 9.6 
= D? — (K^ p°)pr* 7.0 
at pt p )P 
D'2K D? > (Km Ja wavep" 5.1 
D? + (K-ph)vna* 6.8 
_ D? > (Kr )pr 8.5 
D? + K*a*a D? + K”? (s * x )s 8.9 
D? 3 (K*?n )yn* 9.7 
DE => (UC m )S=waveT JAT > 10.0 
D? + K ((1*n )sr*)A > 10.0 
D? — (K^ r*)s-wave(n tr" )s » 10.0 
z = D [S] > (K~1*)v(a*17) 8.8 
D K tt V V 
cree UIS D^ + (K^ jsa t t )y 5.8 
D? — (Kx *)v(n^m )s > 10.0 
D? — (K^ n*)v(n^m )s 6.8 
D? > (K^ )s-wave( ^) 9.7 


VI. SYSTEMATIC UNCERTAINTIES 


'The source of systematic uncertainties are divided into 
four categories: (I) amplitude model, (IT) background es- 
timation, (III) experimental effects and (IV) fitter per- 
formance. The systematic uncertainties of the free pa- 
rameters in the fit and the fit fractions due to different 
contributions are given in units of the statistical standard 
deviations cstat in Tables[VII]IX] These uncertainties are 
added in quadrature, as they are uncorrelated, to obtain 
the total systematic uncertainties. 


1. Amplitude model 


Three sources are considered for the systematic un- 
certainty due to the amplitude model: the masses and 
widths of the K7 (1270) and the aï (1260), the barrier 
effective radius R and the fixed parameters in the Kr S- 
wave model. The uncertainty associated with the mass 
and width of K7 (1270) and the aj (1260) are estimated 
by varying the corresponding masses and widths with lo 
of errors quoted in PDG [1], respectively. The uncer- 
tainty related to the barrier effective radius R is esti- 
mated by varying R within 1.5 — 4.5 GeV! for the in- 
termediate resonances and 3.0 — 7.0 GeV! for the D? 
in the fit. The uncertainty from the input parameters of 
the Kr S-wave model are evaluated by varying the input 


201711.02283v1 


chinaXiv 


10 


TABLE V. Fit fractions for different components. The first and second uncertainties are statistical and systematic, respectively. 


Component Fit fraction (76) Mark III's result E691's result 
D" > Kp 12.3 + 0.4 +0.5 14.2 + 1.6 + 13 Z 2 


D? — K- aï (1260)(p°x*) 


54.6 +2.8+3.7 49.2+24+8 47+5+10 


D? — Kj (1270)(K*"r-)m+ 0.8 + 0.2 + 0.2 


D" > Ky Q270(K p )* 3.4+0.3+0.5 


6.6 + 1.9 +3 - 


D? > K ntp 8.4 + 1.1 + 2.5 


8.4 + 2.2 +4 5 + 322 


D? = Kort ag 7.0 + 0.4 + 0.5 


14.0 + 1.8 + 4 I13E23:E3 


D? 3 K n'ant” 


21.93 0.6 + 0.6 


24.2 + 2.5 + 6 23+2+3 
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FIG. 3. Distribution of (a) TU ts (b) Dus (c) Matr- 


(h) m 


K-mt 


= E 1.6 0 E 
m(K*2*) (GeV/c?) 


; (d) Duis (e) Mk-rfr-’ (£) Mnt? (g) Tt aia and 


+„+, where the dots with error are data, and curves are for the fit projections. The small red histograms in each 
1^2 


projection shows the D? > KO «* peaking background. In (d), a peak of K2 can be seen, which is consistent with the MC 
expectation. The dip around the K9 peak is caused by the requirements used to suppress the D? > K9K-* background. 
Plot (i) shows the fit (curve) to the distribution of the x (points with error bars) with a Gaussian function and the fitted values 


of the parameters (mean and width of Gaussian). 


values within their uncertainties. All the change of the 
results with respect to the nominal one are taken as the 
systematic uncertainties. 


2. Background estimation 


The sources of systematic uncertainty related to the 
background include the amplitude and shape of the back- 
ground D? > K&K~7*, and the other potential back- 
grounds. The uncertainties related to the background 
D? + K9K-7* is estimated by varying the number of 
background events within lo of uncertainties and chang- 


ing the shape according to the uncertainties in PDF pa- 
rameters from CLEO . The uncertainty due to the 
the other potential background is estimated by including 
the corresponding background (estimated from generic 
MC sample) in the fit. 


3. Experimental effects 


The uncertainty related to the experimental effects in- 
cludes two separate components: the acceptance differ- 
ence between MC simulations and data caused by track- 
ing and PID efficiencies, and the detector resolution. To 
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TABLE VI. Phases and fit fractions for different amplitudes. The first and second uncertainties are statistical and systematic, 


respectively. 
Amplitude Qi Fit fraction (%) 
DS] > K*p? 2.35 + 0.06 + 0.18 6.5 + 0.5 + 0.8 
D9[P] > K* p? 2.25 + 0.08 +0.15 2.3 +0.2 + 0.1 
D°[D] > K* p? 2.49 + 0.06 +0.11 7.9 + 0.4 + 0.7 
D? — K~a} (1260), aj (1260)[S] > pr + O(fixed) 53.2 2.8 4.0 
D? — K~a} (1260), aj (1260)[D] > p*r + 2.11 +0.15 +0.21 0.33 0.1 + 0.1 
D? — Kj (1270)r+, Kj (1270)[S] > K*?x- 1.48 +0.21 +0.24 0.10.1 0.1 
D? + Kj (1270)x+, Ky (1270)[D] > K*°r 3.00+0.09+0.15 0.7 0.20.2 
D? > Kj (1270)r+, Kj (1270) > K- p? 2.46 + 0.06 +0.21 3.4 0.3 + 0.5 


D? > (gh K )A 
D? — (K- p°)pa 
D? > K - n )s wave p° 
D9 => K-po)yn 
D? — (K*°n-)pat 

D? + K*%(nta-)gs 

D? — (K*°n-)yat 

D? + ((K-nt)S_wavem Jart 
D? + K-((*n-)sn*)A 

D? + (K-m*)s wave (7*1 )s 


D?[S] > (€^n*)v(n*n-)v 
D°’ > (K T *)s- wave ln a)y 
D? +> (Kort )y(ntr")s 
D? + (K-rt)r(ntr)s 
D? — (K-m*)s wave (7 ta )m 


TABLE VII. Systematic uncertainties on masses and widths 
of intermediate resonances K*° and p°. 


Source (estat) 

I IH II IV 
2.21 0.04 0.13 0.10 2.22 
Ig«o 0.87 0.05 0.17 0.07 0.89 
mM. „o 2.37 0.08 0.12 0.08 2.37 


p 


Loo 1.16 0.04 0.11 0.12 1.17 


Parameter total (Ostat) 


M g«o 


determine the systematic uncertainty due to tracking and 
PID efficiencies, we alter the fit by shifting the ye(p) in 
Eq. (9) within its uncertainty, and the changes of the 
nominal results is taken as the systematic uncertainty. 
The uncertainty caused by resolution is determined as 
the difference between the pull distribution results ob- 
tained from simulated data using generated and fitted 
four-momenta, as described in Sec. 


4. Fitter performance 


The uncertainty from the fit process is evaluated by 
studying toy MC samples. An ensemble of 250 sets of 
SIGNAL MC samples with a size equal to the data sam- 
ple are generated according to the nominal results in 
this analysis. The SIGNAL MC samples are fed into 
the event selection, and the same amplitude analysis is 
performed on each simulated sample. The pull variables, 


Tt, (P K-)A[D] > K=p 


0.43 + 0.09 + 0.12 1.1 + 0.2 + 0.3 
0.14 + 0.11 + 0.10 7.4 1.6 5.7 
2.45 + 0.19 + 0.47 2.00.7 1.9 
1.34 + 0.12 + 0.09 0.4 3: 0.1 0.1 
2.09 + 0.12 + 0.22 2.4 + 0.5 + 0.5 
0.17 + 0.113 0.12 2.60.6 £0.6 
2.13 £0.10 0.11 0.8 0.1 X 0.1 
1.36 + 0.08 30.37. 5.6 3: 0.9 + 2.7 
2.23 + 0.08 + 0.22 13.1 + 1.9 + 2.2 
1.40 + 0.04 + 0.22 16.3 + 0.5 + 0.6 
1.59 +0.13 +0.41 5.4 +1.2+1.9 
0.16 + 0.17 +0.43 1.9 + 0.6 + 1.2 
2.58 É 0.08 E 0.25 2.9 30.5 É 1.7 
2:92 E 0.14 £0.12 0.3 30.1 £0.1 
2.45 + 0.12 + 0.37 0.5 + 0.1 + 0.1 


Vinput — Vii 
P. are defined to evaluate the corresponding un- 


certainty, where Vinput iS the input value in the generator, 
Vat and og, are the output value and the corresponding 
statistical uncertainty, respectively. The distribution of 
pull values for the 250 sets of sample are expected to be 
a normal Gaussian distribution, and any shift on mean 
and widths indicate the bias on the fit values and its 
statistical uncertainty, respectively. 


Small biases for some fitted parameters and fit frac- 
tions are observed. For the pull mean, the largest bias is 
about 1996 of a statistical uncertainty with a deviation 
of about 3.00 from zero. For the pull width, the largest 
shift is 0.87 + 0.04, about 3.0 standard deviations from 
1.0. We add in quadrature the mean and the mean error 
in the pull and multiply this number with the statistical 
error to get the systematic error. The fit results are given 
in Tables[XHXTI] The uncertainties in Tables[X] XII] are 
the statistical uncertainties of the fits to the pull distri- 
butions. 


VII. CONCLUSION 


An amplitude analysis of the decay D? > K-n*ntn 
has been performed with the 2.93 fb ^! of e*e- collision 
data at the (3770) resonance collected by the BESIII de- 
tector. The dominant components, D? + K~af (1260), 
D? — K*°p°, D? => four-body nonresonant decay and 
three-body tnn D? — K-«*g? improve upon 
the earlier results from Mark III and are consistent with 
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TABLE VIII. Systematic uncertainties on fit fractions for different components. 


Source (Ostat) 


Fit fraction I IH II IV 


total (Ostat) 


DU > K*0p0 1.12 0.06 0.11 0.08 1.13 
D? — K-at (1260) 1.32 0.09 0.12 0.06 1.33 


D? — Kj (1270) (K*9z—)n* 1.41 0.02 0.12 0.10 1.42 
D? — Kj (1270)(K~ p°)x+ 1.58 0.04 0.23 0.06 1.60 


D? — K~ xt p? 2.22 0.10 0.12 0.15 2.23 
D? + K*og*z- 1.32 0.08 0.13 0.10 1.34 
D? > K-ntata- 0.94 0.10 0.09 0.12 1.00 


TABLE IX. Systematic uncertainties on phases and fit fractions for different amplitudes. 


Source (estat) 


$i I II mp gv total (Ostat) 
D'S] > KF 2.96 0.04 0.14 0.13 2.97 
D? [P] > K*9„0 1.98 0.04 0.11 0.12 1.98 
D9|D] > K*090 1.78 0.03 0.18 0.09 1.79 


D? — K-a} (1260), a} (1260)[D] > p r+ 1.38 0.02 0.09 0.09 1.39 
D? — Ki (1270)&* ,K, (1270)[S] > K*°r~ 1.10 0.07 0.10 0.09 1.11 
D? — Ky (1270)&*,K, (1270)[D] > K*°r- 1.61 0.06 0.11 0.06 1.62 
D? ^ Ky (1270)a+,Ky (1270) > K~p® 3.61 0.03 0.09 0.13 3.62 


D’ > (WK -)ant 1.28 0.06 0.14 0.09 1.29 
D? — (K~p°)pat 0.92 0.10 0.10 0.07 0.93 
D? — (Krt )s—-wavep? 2.46 0.06 0.10 0.09 2.47 
D? — (K7 pP )yrt 0.74 0.01 0.09 0.08 0.75 
D? — (K*95—)pn* 1.82 0.03 0.09 0.06 1.82 
D? + K*(n*v-)s 1.07 0.04 0.12 0.11 1.08 
D? — (K*9g—)ymt 1.00 0.02 0.10 0.18 1.02 
D? — ((K—m*)s wave )AT T 4.78 0.15 0.12 0.07 4.79 
D? — K-((n*n-)sgn*)A 2.69 0.13 0.10 0.07 2.70 
D? => (K—m*)s wave (1 ta )g 6.27 0.04 0.10 0.12 6.27 
D°[S] > (K-nt)y(nta-)y 3.28 0.06 0.09 0.06 3.28 
D? — (K~2+)g_wave(mta-)v 2.59 0.09 0.10 0.10 2.60 
D? + (K-tt)v(nmta-)s 3.07 0.09 0.10 0.18 3.08 
D? + (K-n*)qv(n**-)s 0.81 0.04 0.12 0.06 0.82 
D? — (K—m*)s wave(at 1^ )m 3.11 0.06 0.11 0.16 3.19 


Source (estat) 


1 g mm Iy fetal (Ostat) 


Fit fraction 


DIS] > Kp 1.76 0.04 0.09 0.10 1.77 
D^ [P] > K*959 0.27 0.02 0.09 0.12 0.31 
D°[D] > K*950 1.79 0.06 0.12 0.17 1.80 


D? + K-aj (1260), aï (1260)[S] > pz — 1.48 0.10 0.12 0.07 1.45 
D? — K-aï (1260), aï (1260)LD] > ox 0.93 0.04 0.09 0.06 0.94 
D? — Ky (1270)&*,K, (1270)[S] > K*?x-— 1.01 0.05 0.11 0.16 1.03 
D? > Kj (1270)n* ,K, (1270)[D] > K*%r- 1.12 0.03 0.12 0.13 1.14 
D? — K4(1270)- 4 *,K, (1270) > K - p? 1.58 0.04 0.23 0.06 1.60 


D? — (p K )Am t 1.38 0.08 0.09 0.09 1.39 
D? — (K*?n)pm 0.93 0.06 0.09 0.16 0.95 
DE — (Kn )s wave p? 2.81 0.09 0.11 0.09 — 2.82 
D? => (K-ph)yn* 0.69 0.03 0.09 0.06 0.70 
D? — (K*°r-)pat 0.93 0.06 0.09 0.16 0.95 
D? + K*9 (ntr )s 1.06 0.05 0.09 0.20 1.08 
D? — (K*9n—)yn* 0.60 0.02 0.00 0.10 0.61 
D? — ((K—n*)s wave Jant 3.10 0.07 0.09 0.06 3.10 
D? => K-((n*n-)sn*)A 1.14 0.08 0.10 0.07 1.15 
D? — (Km )s wave * 7 )s 1.29 0.12 0.10 0.12 1.30 
D9?[S] > (K-c*)y(r*n-)v 1.73 0.07 0.09 0.07 1.73 
D? — (K—m*)s wave(1^1-)v 2.08 0.12 0.10 0.07 2.09 
D? => (K-m*)vy(n*7-—)s 3.54 0.05 0.10 0.11 3.54 
D? = (K-n*)q(n*7-)s 0.87 0.07 0.11 0.07 0.88 
D? — (K—m*)s wave (1^1) 0.99 0.09 0.10 0.08 1.01 
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TABLE X. Pull mean and pull width of the pull distributions for the fitted masses and widths of intermediate resonances K*? 
and o? from simulated data using either the generated or fitted four-momenta. 
Generated pi Fitted p; 
pull mean pull width pull mean pull width 
^ mgo  0.07+0.07 1.05 30.05 0.06 +0.07 1.04 0.05 
Ig«o —0.03 + 0.06 0.97 + 0.04 —0.17 + 0.06 0.97 + 0.04 


Parameter 


mo 0.03 + 0.07 1.06 +0.05 —0.02 + 0.07 1.06 + 0.05 
Io 0.10 +0.07 1.08 +0.05 0.06 3- 0.07 1.07 + 0.05 


TABLE XI. Pull mean and pull width of the pull distributions for the different components from simulated data using either 
the generated or fitted four-momenta. 


: 3 Generated p; Fitted p; 
Fit traction pull mean pull width pull mean pull width 
DU > K*050 0.05 +0.06 0.92+0.04 0.04 3- 0.06 |. 0.89 4- 0.04 
D’ > K- aï (1260) 0.02 +0.06 0.91 +0.04 0.04 + 0.06 0.87 + 0.04 


D? — K7 (1270)(K*°n-)r+ —0.08 + 0.06 0.98 +0.04 —0.06 + 0.06 0.97 + 0.04 
D? — Kj (1270)(K-p")r+  0.0012-0.06 0.98+0.04 0.01-- 0.06 0.99 0.04 


D? — K~xt p 0.142- 0.06. 0.92 -- 0.04. 0.11 20.06. 0.88 4- 0.04 
D? => K*Ogta- —0.08 + 0.06 0.96 + 0.04 —0.09 + 0.06 0.96 + 0.04 
D? => K-ntata- 0.10+0.06 0.9440.04 0.12+0.06 0.93+0.04 


TABLE XII. Pull mean and pull width of the pull distributions for the phases and fit fractions of different amplitudes, from 
simulated data using either the generated or fitted four-momenta. 


Generated pi Fitted p; 
$i pull mean pull width pull mean pull width 
DS] > KP 0.11 £0.06 1.01£0.05 0.08£0.06 1.00 0.04 
D°[P] > K*990 0.10 + 0.07 1.03 2- 0.05 0.08+0.06 1.02 4 0.05 
D°[D] > K*959 0.05 +0.07 1.04 -- 0.05 0.01--0.07 1.03 0.05 


D? ^ K- af (1260), aï (1260)[D] > p'z* —0.07 + 0.06 1.02 + 0.05 —0.05 + 0.06 1.02 + 0.05 
D? > Kr (1270)&*,K, (1270)[S] > K*°,- 0.06 0.07 1.03 +0.05 0.06 0.06 1.03 4 0.05 
D? > Kj (1270)r a m E D = —0.02 + 0.06 0.98 + 0.04 —0.06 + 0.06 0.97 + 0.04 
D? +> xT an F Kz (1270) > 0.12+0.06 1.00+0.04 0.11+0.06 1.00 0.04 


D? > —0.06 + 0.07 1.05 + 0.05 —0.09 + 0.07 1.05 + 0.05 
D? + K- p É + —0.03 + 0.06 0.96 + 0.04 —0.01 + 0.06 0.96 + 0.04 
Do =$ K= H) s emit? - 0.07 + 0.06 0.92 + 0.04 —0.08 + 0.06 0.92 + 0.04 
D? — (K~p°)yat —0.05 + 0.06 1.02 + 0.05 —0.07 + 0.06 1.01 + 0.05 
D? — (K*9n pat 0.00 +0.06 0.99 -- 0.04 0.00 «c 0.06 0.99 + 0.04 
D? + K*°(ata7)s —0.08 + 0.07 1.03 +0.05 —0.11 + 0.07 1.03 + 0.05 
D? — (K*9g yat 0.17 +0.06 0.99 --0.04 0.15 +0.06 0.98 + 0.04 
D? — ((K-n*)s wave )AT —0.04 + 0.06 0.92 + 0.04 0.02 4 0.06 0.92 + 0.04 
D? — K- ((1 *n-)sn*)A 0.00 + 0.07 1.05 + 0.05 —0.02 + 0.07 1.04 + 0.05 
DY => Kt) wave( t r" )s 0.10 + 0.06 0.98 +0.04 0.08+0.06 0.98 + 0.04 
D? [S] > (K-«^)v(r**-)v —0.02 + 0.06 0.97 + 0.04 —0.03 + 0.06 0.98 + 0.04 
D? — (K-tt)g-_wave(n ta )v 0.08 + 0.06 0.93+0.04 0.06 +0.06 0.92 0.04 
D? — (K-m*)y(1*7—)s —0.17 + 0.06 0.94 + 0.04 —0.17 + 0.06 0.94 + 0.04 
D? — (K-tt)T(nmta-)s 0.01 + 0.06 1.01 + 0.05 —0.02 + 0.06 1.00 + 0.04 
D? — (K-rt)s-wave(nt r )r 0.14 +0.07 1.12+0.05 0.12+0.07 1.11 0.05 
" " Generated p; Fitted p; 
Pee pull mean pull width pull mean pull width 
D'[S|] > K*V pU 0.08 £0.06 0.88 £0.04 0.07 £0.06 0.87 £0.04 
D°[P] + K*°p° 0.10+0.06 0.97 +0.04 0.10+0.06 0.96 + 0.04 
D°[D] => K*°p° —0.15 + 0.07 1.10 +0.05 —0.15 + 0.07 1.10 + 0.05 


D? — K~a} (1260), aj (1260)[S] > prt 0.03 +0.06 0.91 +0.04 0.04 +0.06 0.90 0.04 
D? — K~a} (1260), af (1260)[D] > p°x+ 0.02+0.06 1.00 +0.04 0.03 +0.06 1.00 0.04 
D? — Kj (1270)n,K] (1270)[S] > K*9x- —0.14 + 0.07 1.02 + 0.05 —0.18 + 0.07 1.09 + 0.05 
D => Kj (1270)&*,K, (1270)[D] + K*9x- —0.11 +0.06 0.99 + 0.04 —0.09 + 0.06 0.99 + 0.04 
D? — K; (1270)&^,K, (1270) > K- p° 0.01 + 0.06 0.98 2-0.04 0.01+0.06 0.98 + 0.04 


D? — (p K-)amt 0.06 +0.06 1.00 +0.04 0.04 3- 0.06 0.99 0.04 
D? — (K7 p )prt 0.11 +0.06 0.95 +0.04 0.0940.06 0.94 0.04 
D9 — (Km )s wave p? 0.05 +0.07 1.04 3- 0.05 0.05 3- 0.07 1.04 0.05 
D? > (K-pt)yn + 0.01 +0.06 0.98 +0.04 0.02 3- 0.06 0.97 0.04 
D? — (K*9n )pn* 0.15 +0.06 0.93 +0.04 0.15 +0.06 0.93 + 0.04 
D? + K*9 (zt x—)s —0.19 + 0.06 1.03 + 0.05 —0.18 + 0.06 1.02 + 0.05 
D? + (K*9g-)ymt —0.08 + 0.06 1.00 + 0.04 —0.09 + 0.06 1.00 + 0.04 
D? — ((K-nt)8_wavemr )AT T 0.02 +0.06 0.98 +0.04 0.02 4- 0.06 0.97 + 0.04 
D? > K- ((n*2-)sn*)A 0.04 + 0.06 1.01 +0.05 0.04 +0.06 1.00 4 0.04 
D? > (K-m*)s wave (7 * r" )s —0.10 + 0.06 0.93 + 0.04 —0.09 + 0.06 0.93 + 0.04 
D°[S] > (K-c*)y(n*2-)vy 0.03 3-0.06 1.02 +0.05 0.03 +0.06 1.01 0.05 
D? + (K-m*)s wave (7 f )y 0.04 + 0.06 1.00 +0.04 0.04 3- 0.06 0.99 0.04 
D9 > n T)v(*7-)s 0.09 3-0.07 1.06 +0.05 0.1140.07 1.04 4 0.05 


Km) 
D? 5 (Kn )p(n*2-)s 0.01 +0.07 1.0540.05 0.00 +0.07 1.03 + 0.05 
D? — (K-nt)g_wave(nta-)r 0.05 +0.06 0.96 +0.04 0.05 +0.06 0.96 + 0.04 
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them within corresponding uncertainties. The resonance 
Ky (1270) observed by Mark III is also confirmed in this 
analysis. The detailed results are listed in Table [V] 


About 4096 of components comes from the nonreso- 
nant four-body (D° — K-m«*m«*-«- ) and three-body 
(DU > Korto and D? => K*°rtn—) decays. A de- 
tailed study considering the different orbital angular mo- 
mentum is performed, which was not included in the 
analyses of Mark III and E691. An especially interest- 
ing process involving the Kr S-wave is described by an 
effective range parametrization. 


By using the inclusive branching fraction B(D? — 
K-xtntn-) = (8.07 + 0.23)% taken from the PDG 
and the fit fraction for the different components F'F(n) 
obtained in this analysis, we calculate the exclusive ab- 
solute branching fractions for the individual components 
with B(n) = B(D? > K-atntn-) x FF(n). The re- 
sults are summarized in Table and are compared 
with the values quoted in PDG. Our results have much 
improved precision; they may shed light in a theoret- 
ical calculation. The knowledge of D? > K*°p° and 
D? — K~a} (1260) increase our understanding of the de- 
cay D? — VV and D — AP, both of which are lacking 
in experimental measurements, but have large contribu- 
tions to the D? decays. Furthermore, knowledge of the 
submodes in the decay D? > K^z*«*«- will improve 
the determination of the reconstruction efficiency when 
this mode is used to tag the D? as part of other mea- 
surements, like measurements of branching fractions, the 
strong phase or the angle 7. 


TABLE XIII. Absolute branching fractions of the seven com- 
ponents and the corresponding values in the PDG. Here, we 
denote K*? > K-nt and p? > «*«-. The first two un- 
certainties are statistical and systematic, respectively. The 
third uncertainties are propagated from the uncertainty of 
B(D9 > K-atatn). 
Component 


Branching fraction (7%) PDG value (96) 


DV =>. KY 0.99 + 0.04 + 0.04 + 0.03 1.05 + 0.23 
D? + K-alt(1260)(px*) 4.41 +0.22 +0.30+0.13 — 3.6 + 0.6 
D > K] (270) nat 0.07 + 0.01 + 0.02 + 0.00 955 + 0.03 
D? > KT (1270)( K7 p’)x* 0.27 + 0.02 + 0.04 + 0.01 E 
D? > K-«*p? 0.68 + 0.09 + 0.20 + 0.02 0.51 + 0.23 
D? + K*ogta- 0.57 + 0.03 + 0.04 + 0.02 0.99 + 0.23 
D? + K-atatn- 1.77 + 0.05 + 0.04 + 0.05 1.88 + 0.26 
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VIII. APPENDIX A: AMPLITUDES TESTED 


The amplitudes listed below are tested when deter- 
mining the nominal fit model, but not used in our final 
fit result. 

Cascade dk 


Ne K-)m i ac D-wave 
AVIT —)r*, K*?n- S and D-waves 
*- (1410)(K*92— m * 
27 (1430)(K*95—)n*, K2- (1430)(K - p)n* 
*- (1680) (Ks), K*- (1680) (K - ph)n* 
2 (1770)(K*'n ~ K3 (1770)(K~ p°)xt 


ios 30) 

Kn (1300)(p°x*) 

K-at (1260)(fo(500)x*) 
Quasi-two-body amplitudes 
K*° fo (500) 

K*° fo(980) 

Three-body amplitudes 
K*9(n*n-)y S, P- and D-waves 
(K-n*)vp? S, P and D-waves 
K3 (1430) (7 —)s 
K39(1430) p° 

K” fo(1270) 

(Kv )s fo(1270) 


K*95—-)Am*, K*?n- S and D-waves 
Four-body nonresonance amplitudes 
(K^ m*)r(n*n-)vy P- and D-waves 


201711.02283V1 


chinaXiv 


7 )r P- and D-waves 


7 )y P- and D-waves 
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